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TECHNICAL PAPER 


mode-medium instability and its correction with 

A GAUSSIAN REFLECTIVITY MIRROR 


I. INTRODUCTION 


The mode-medium instability (MMI, effect is due to 
the upper and lower CO, fasttrig levels which 'bating rate in an unstable laser resonator with hard- 

where the resonator intensity is g * • wholly to the radiation intensity variations pro- 
edge output mirror is attributed mainly, although not ^ rate leads to acousti c density 

duced by diffraction from the edge of the output mirro . dielectric constant which in turn causes the 

perturbations and, thus, spatial and temporal ^^^ variati<ms in the heatin g rate 

radiation to be diffracted even more. The change t _ y degrades with time, and 

and the feedback cycle is complete^ As ^ be co-p^ely destroyed. Thus, MMI 

for pulses of sufficient intensity and duration th q y n . A rev i ew of the earlier experi- 

is of major importance in designing high-power long-pulse CO, lasers. A review 

mental work is given in References 1 through 8. 

The numerical modeling of MMI , evolves fhe soluuon offwo “ 

rial equaiions as shown later. Recently, in the „Ltable resonator. The computation, 

used to study the qualitative time evolution exDected larger Fresnel number and smaller 

have prompted many suggestions lor improvement. 

The suggestions for unproving beam qual.ty in** 

[6], pressure reduction 16], and phase conjugate * ion rj pp les pro duced by the edge of the output 

mirror designs aimed at reducing or e lmina ing characterized by a step phase discontinuity which 

mirror of the unstable resonator. The steppe m . interfere with and cancel the other set of 

produces an additional set of edge waves that ar - g ^ other hand the so ft-edge mirror is 

edge waves produced by the outer edge of t e . diffraction produced by the hard-edge 

characterized by a tapered edge reflectivity which mirror due to its 

output mirror. The Gaussian reflectivity mirror [13-17 « into a smooth Gaussian 
ability to transform the diffraction modu a e mo „ , output bea m of the unstable resonator. A 

mode without sacrificing the large mo e ’ i° however is not available. Another suggestion, known as 
quantitative investigation, including MM \\ h ’ ‘ the 5eam aroU nd. This type of possible 

beam dithering, is to average out the heat. ^ _ * beam must be moved on the order of micro- 

solution would be very difficult to implementsince the beam quality. It has been 

second scales. Pressure reduction is also an imp of Lr-f.eld energy delivered to 

demonstrated [6] that a reduction m pressu expense of total output power. For higher 

the far-field. However, as will be shown later, this is at the expense 



power it is still necessary to operate at higher pressure. It has also been suggested that phase coniueation 

O nnotiva't I'" ‘° CanCe ' medi,,m indUCed Phase iis '° nims ’ to reduce the iSXT 

9 »n T T 1 g ' Ve 3 qlU ' nm: " lvc slud y of MMI >™ small Fresnel number to complement Reference 
9 and to study the corrective mechanisms suggested previously by others. Kelerence 

HcX Th COmP d mm0T a " d an “ CtiVe "° nlinear medium causes a great deal of numerical d f 
w"h Lire .Sn sS‘ !n U ,heTo„ a eh 7*7f* e ‘ ,Ua,i ° nS " " l " d h <* dimension 1"' 

’ZZZ*- The use of a ? tional 8ain sheets improws -* SKE 

^KSSSKSbs 

t( wfth T Th.s paper reports the first time resolved mode calculations in an unstable resona 

w.th a Gaussian reflectivity mirror and with MMI effects included in the calculation. 

n ln t lhe followin 8 sections, the standard formalism, for both the hard-edge and Gaussian 
g,vt y i„ "heT„d POWer ° U ' Pm are Prese "' ed in ,he f0Urth A discussio " concluding the paper Z 


II. GENERAL FORMALISM 


In the active medium, the dielectric constant can be written as [9] 


— 1 + Ae = I + 4 Tui + jo -11 

nh) * OJ 


( 1 ) 




8 = 


go 


1 +. 


( 2 ) 


derived* etnlier hT" Si8na ' 8ai " ““ Sa,Ura,i ° n inKnSity ' A " eqUa,ion of s,ate for P ^ been 
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Figure 1. Three-level lasing system. 
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where C, is the speed of sound, y = 1 .7 is the ratio of specific heat at constant pressure to the specific heat 
at constant volume, and q is the heating rate per unit volume which is related to the intensity / by [9] 


V’l 


gl + \Y) ^ xal ~ 8 ^ x ‘ 


(?) 




i+l»l 


(4) 


where equation (2) has been included and / has been normalized according to 


\u\ 2 = 


(5) 


The solution of equation (3) is [4] 


-('V — 1 ) r 1 

&*■» = T x J 

q(x,t')-}r q(x + CJt — t'),t' ) — }- q{x — C s (t t ),t ) 

c ; o 

2 2 


Equations (1) through (6) are to be solved with the Maxwell equation for the field amplitude, 


E = E exp — Ac)] 


(7) 


where E satisfies the envelope equation. 
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The wave number A. and the frequency oo are related by oj = ck, and 


/ = 



(9) 


The set of equations given in (I) to (9) would be a satisfactory model for the MM1. In actual 
execution of the solution, the following approximations and limitations are introduced: (i) The model is 
limited to one transverse dimension which does not change the qualitative significance of the results, 
(ii) The Fast Fourier Transform with 256 data points is used. This limits the width of the system and the 
spatial frequency content of the transverse dimension as shown later, (iii) The gain and density perturba- 
tions to the dielectric constant in (1) are applied as complex gain sheets at discrete intervals up and down 
the resonator cavity. This approximation may only affect the quantitative accuracy of the results, 
(iv) Equation (6) is simplified by expanding q in a Taylor series about jc and using a summation to 
approximate the time integration. These approximations are numerically accurate as long as the distance 
between data points and the time increment in the summation are small. 


III. NUMERICAL PROCEDURES 


A. Free-Space Propagation 

In this work, a Fast Fourier Transform (FFT) [20] technique is used to repeatedly solve the free- 
space wave equation up and down the resonator cavity. The medium effects are applied as complex gain 
sheets at discrete intervals throughout the cavity. For simplicity, the calculation is limited to one 
transverse dimension. All physical constants are given in their normal three-dimensional units. The total 
output power is given in units of power per unit length in the other dimension. Any explicit time 
dependencies are handled by replacing t with 


t 


tn = 



( 10 ) 


where n is the number of trips around the resonator. To begin, equation (9) is rewritten in a dimensionless 
form for a fixed time interval, t n . 




Av 2 dz 


u(x,z,tfj) = - 2nxiNfi c (x,z,t n )u(x,z,t„) 


( 11 ) 
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where x is scaled by the transverse radius, a, of the output mirror, and z is scaled by the length, L , ot the 
cavity. The Fresnel number, Nj, is defined by 


Nr = 


a 

LK 


(12) 


and the envelope, E, of the electric field is normalized according to 


I P 

\ur = 




(13) 


Here the dielectric constant of (1), which contains all of the medium effects, is rewritten as a complex 
gain profile. 


G c {x,z,t„) = -ikL&e(x,z,tn) 


gL — 4i:ikLa 



(14) 


where the real part is the gain and the imaginary part is the phase shift due to the density variation. 

In order to solve (11), the standard gain sheet approximation [18] is assumed by dividing the 
cavity into N s sections, as depicted in Figure 2, and free-space propagation through each of these sections 
is assumed. All of the medium effects accumulated by the beam as it passes through a section are lumped 
together into a single, transverse complex gain profile (sheet), which is applied to the beam after tree- 
space propagation through a section. To accomplish this approximation, G r is assumed to vary slowly 
along the axis of the resonator by making the following substitution: 


u(x,z,tn) = U(x,z,t n ) exp 


m 


s! 


I 

a 




N 


M‘ 


Z4 z 3 Z 2 Z| Zq 


Figure 2. Resonator divided into N s sections. 


(15) 
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Equation (II) now takes the form of a free-space wave equation. 


-^-4 iriNj-S- 
Ar 1 


dz I 


U(x,z,tff) = 0 


( 16 ) 


Taking the FFT, equation (17) is obtained 


c3(7(k,z,^) 

dz 


— <K~ 

4ixN r 


U(K,Z,t n ) 


(17) 


where 


0(k, z,t Tl ) = FFT \U{x,z,tn) 1 


(18) 


Equation (17) has the solution 


U(K,z,t n ) = U{K,z 0 ,tn) exp 




(19) 


Thus, the propagation of a wave front through a section is given by 


U(x,z m ,tfj) = FF'T 1 


exp(^ 

F \4ttN, 


- Az) 


FFT f t/(.r , z„, i , /, 


,T)ij , 


( 20 ) 


where z m ^ and z m are the longitudinal positions of the front and back edges of the m " section, and 


Az = z m — z, 


m *m— I 


A h 


( 21 ) 


Remember that the z-dimension has been normalized to the length of the cavity. 
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B. Complex Gain Profile 


After free-space propagation through a section, the medium effects are applied as a complex gain 
sheet which, according to (15), yields 


u(x,z m ,tjt) = exp ^ i G c (x,z n „tfi)kzj x FFT 1 

The use of complex gain sheets is valid as long as G c (x,z„„t„) varies slowly over the distance through a 
section. Since G c contains the medium effects, its size should match that of the lasing medium inside the 
laser cavity. In our calculations, it is assumed that the medium fills the entire cavity and G< covers the 
entire calculation space. The time dependence of G< is accommodated by using the field protile, u(t^-\), 
of the previous trip to calculate a new complex gain profile, G ( .(7, T ), for the next trip. In essence. 




Now G ( . is broken down into its real and imaginary parts. 


eXP (^ A 7 


FFT[u(.v,z„,_i,G)] 


] 

_ 


[ 22 ) 


G ( = G, + iGj , 


(24) 


and each part is calculated independently. According to (14) and (2), the real part of the complex gain, 
G,., depends explicitly on the intensity and is given by 


G,(x,z m , tfj ) 


So 

IT U U,Z,„,/,t)|" + U (x,Z„,,tfj)\ 


(25) 


where the + and — superscripts denote the forward- and backward-traveling wave. 

Calculation of the imaginary part, G„ of the complex gain profile is not so straightforward since it 
involves solving the time integral of (6). Equation (6) is substituted into (14) to obtain 


V m 0 C~ s f 0 


q{x,z m ,t ) — c/(.\ + C s (t, t t ) 


- i q(x-C x (t n - t'),z„„t') 


(26) 
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The integral is evaluated by expanding q in a Taylor series about x, substituting a summation for the 
integration, and transforming the summation into a form that can be computed with a running total, rather 
than performing the summation every time. The second derivative of q with respect to x resulting from 
the Taylor series expansion is computed with finite-difference differentiation formulas [21]. As 
mentioned before, the accuracy of these approximations is good as long as the distance between data 
points and the time increment in the summation are small. Appropriate increments are determined by 
repeatedly performing the calculations with smaller and smaller increments, until the desired accuracy is 
obtained. In our calculations, N s = 10 complex gain sheets, 256 transverse data points, and a time 
increment equal to the round-trip propagation time of light are used. 

C. Feedback and Output Mirrors 

After the EM wave has been numerically propagated from one end of the resonator to the other, it 
is necessary to model the reflection of the wave from the end mirror. For the hard-edge unstable resona- 
tor, it is assumed that the end mirrors are perfectly reflecting and spherical. Reflection from a spherical 
mirror induces a phase shift in the wave front. This phase shift is given by [22] 


<]> V /? 2 at] [ 1+ ^2-+ ] - (27) 

where R m is the radius of curvature of the mirror, and <[> has been expanded with a Maclaurin series in x. 
R,„ is taken to be positive for a concave mirror and negative for a convex mirror. In the paraxial limit 
(jr/4 R m <<l), reflection from the feedback mirror (mirror 1) is given by 


u + (x,\,n) = u (x,\,n) (28) 

For the output mirror (mirror 2), a similar expression with only a minor modification to account 
for the output coupling of the laser is used. Because the output is taken outside the edge of the output 
mirror, there is no reflected wave beyond this edge. Therefore, reflection from the output mirror is given 
by 


m + (a,0,«+ I) 


u (a,0,w) expf-^ 
0 


|x | 2 < a 

M 2 > « 


(29) 


where a is the transverse radius of the output mirror. Diffraction is neglected beyond the edge of the 
feedback mirror because it is small in comparison with the output coupling of the laser. 
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Replacing the hard-edge output mirror with a Gaussian-reflectivity mirror (GRM) is accom 
plished by replacing (29) with 


u + {x,0,n + 1) 


= u~(x,0,n) VR(x) exp 



(30) 


where the reflectivity, R( x), of the GRM is given by [16] 


R(x) = Ro exp 



(31) 


where a is the Me of the GRM and R n is the reflectivity at the center {x - 0). 


D. Numerical Problems 

Several numerical problems, including aliasing, energy leakage, and Gibbs' phenomenon, must 
be accounted for when using the FFT to model wavefront propagation. These problems have been 
addressed in a previous work |I8|. For the sake of completeness, aliasing and energy leakage are briefly 
discussed. Gibbs' phenomenon is discussed more extensively because oi an associated problem not notec 

in Reference 18. 

Aliasing occurs when the number of data points is too few for the FFT to adequately r ^ s0 ' v ^ the 
frequency content of the wavefront. For the hard-edge resonator, aliasing limits the use ot the FFT to 
small Fresnel numbers, since the spatial frequency content of the wavefront is proportional to the Fresnel 
number For large Fresnel numbers, the number of data points needed to sample the high frequency con- 
tent of the wavefront becomes prohibitive. By trial and error, it is determined that 256 data points are 
adequate for the small Fresnel number results presented here. 

Next, the energy leakage problem is addressed. Unlike the continuous Fourier transform, whose 
integration limits are infinite, the FFT’s integration limits are finite. Thus, the FFT actually solves the 
problem of a cyclic function, where the sampled function, in this case the wavefront inside the resonator 
represents one complete cycle of the cyclic function. The guard band separating the apertures of each 
cycle must be large enough to prevent diffracted energy from leaking into adjacent cycles as the solution 
is propagated down the resonator. Assuming a uniform plane wave illuminating an aperture of width 2a 
and using Huygens’ integral to estimate the fractional energy falling outside a band of-G^a ~£ * - U b a, 
Sziklas and Siegman [18] developed a criterion for the guard band ratio which is given by 


G h 


I 

* 2tt 2 A r t e b ' 


(32) 
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where e,, is the fractional energy spill-over that can be tolerated, and N c is the collimated Fresnel number 
which is given by 


N,. = M ~ 


M + 



_ Kf 

M 4- I 


n = 2ML n 
M~-\ cq 


(33) 


In theory, equation (32) gives a reasonable estimate for G b . However, in practice, G b must be determined 
by trial and error. If G h is too small, then energy leaking from adjacent cells will produce undesirable 
modulations in the wavefront. Table 1 compares values of G b obtained from (32) with values obtained by 
trial and error as a function ot N e(/ . Here, a fractional energy leakage of e b = 0.004 and a magnification 
of M = 1 .4 are assumed. In practice, G b needs to be larger than the minimum value obtained from (32). 

Finally, Gibbs' phenomenon [23] is addressed. When a discretely sampled wavefront, having a 
sharp cut-off due to a hard-edge aperture, is transformed with the FFT and then back-transformed with 
the inverse FFT, the original wavefront is reproduced exactly. However, if the transformed wavefront is 
propagated a short distance and then back-transformed, a numerical inaccuracy, known as Gibbs’ 
phenomenon, is introduced. When a apertured function is transformed, there will be high frequency 
components due to the sharp edges of the function. Since the number of data points must be finite, the 
FFT effectively truncates the high frequency components at the highest frequency of the transform. In 
essence, the FFT acts as a square low-pass frequency filter. Thus, when the transformed wavefront is 
propagated a short distance, some inaccuracy will be introduced. The iterative application of the FFT, 
propagation, and inverse FFT to perform the mode calculation compounds the problem. 

To reduce or eliminate Gibbs phenomenon, it is necessary to replace the square, low-pass filter 
with a filter that has a more gradual cut-off frequency. There are many different filters that can be used 
lor this purpose [23], This work employs a raised cosine filter [18], also known as a Hanning filter [23], 
depicted in Figure 3(a). By experimenting with several variations of this filter, it was determined that the 
rounded comers of the filter are primarily responsible for the correction of Gibbs’ phenomenon. For the 
hard-edge resonator calculations performed here, it was found that the truncated filter of Figure 3(b) 
yields an accurate solution with the fewest data points. Another filter which may be even better is an 
extended cosine-bell filter [25], where a raised cosine filter is applied to the first and last 10 percent of the 
data and a weight ot unity is applied between. All ot the hard-edge resonator calculations are performed 
with the truncated Hanning filter of Figure 3(b). 

The correction of Gibbs' phenomenon introduces another problem which has not been addressed 
betore. The application of the frequency filter introduces an artificial power loss which is compounded 
by the tact that the tilter is applied several times up and down the cavity; twice for each complex gain 
sheet in the cavity. In a previous work [18], the artificial power loss was not significant since only one 
gain sheet was used. However, in this work 10 gain sheets are used, and the filter is applied 20 times. 
The power loss does not affect the qualitative significance of the results, but a comparison for different 
Fresnel numbers becomes meaningless. When the Fresnel number of the hard-edge resonator is 
increased, the mode volume increases and, as a consequence, the total output power should also increase. 
However, the trequency ot the modulations in the mode increases at the same time, resulting in an 
increased power loss due to the frequency filter. At smaller Fresnel numbers, the increasing mode 
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TABLE 


Figure 3. 


I COMPARISON OF VALUES OF G h CALCULATED FROM EQUATION (32) 
AND OBTAINED BY TRIAL AND ERROR 


N cq 

Calculated 

Trial and Error 

l 

4.0 

4.0 

2 

2.5 

3.0 

3 

2.0 

2.5 

4 

1.8 

2.3 

5 

1.6 

2.0 



Variations of the Hanning filter: (a) normal and (b) truncated. The maximum frequency of 
FFT is proportional to the number of data points. 





volume wins out over the increasing filter losses and the power, output increases. At larger Fresnel 
numbers, the filter losses win out and the power output decreases. Thus, the quantitative inaccuracies are 
greater at larger Fresnel numbers. 

The following is done to correct the artificial power loss: (1) determine the percentage of energy 
removed by the filter, and (2) evenly distribute the loss back into the wavefront. This procedure turns out 
to be important for quantitative comparisons [26], 


IV. RESULTS 


The numerical model developed in the previous section is used to model the time evolution of the 
mode ol the unstable resonator; lirst with a conventional hard-edge output mirror, and then with a correc- 
tive Gaussian-edge output mirror. Figure 4 presents the time history of the intensity and phase profile just 
before the output mirror lor the hard-edge resonator. At 2 |xs the heating has very little effect on the mode 
of the resonator, and the plots of Figure 4(a) are representative, in shape, of the mode of the equivalent 
passive resonator. As the heating builds, the number of diffraction ripples in the intensity profile double, 
and the intensity becomes more strongly peaked near the axis of the resonator compared with the 
intensity at the edge ot the output mirror (.v= 1 ). In addition, the peak-to-peak swing of the phase profile 
grows from about 40° at 2 p,s to about 145° at 1 8 p.s. After 20 |xs the mode becomes very unstable and no 
useful cavity oscillation is obtained. 

A quantitative comparison is begun by varying the equivalent Fresnel number, N ( , q [24], which is 
an important parameter in determining the mode of a hard-edge unstable resonator and is roughly equal to 
the number ol modulations in half of the intensity profile. The dependence of the MMI on N CI/ is 
somewhat complicated because the structure of the resonator mode is determined by diffraction from the 
edge ot the output mirror. As N C(/ is increased, the number of modulations in the mode increase. Thus, 
the MMI ellects are expected to increase as well. However, there is also a corresponding decrease in the 
amplitude ol the modulations which may compensate for the increased number of modulations. Another 
factor affecting the dependence of the MMI on N e(f is related to the amount of power extracted from the 
resonator. Due to the highly diffractive nature of the unstable resonator, the round-trip losses of any 
given mode do not vary consistently with increasing values of N et/ . It has been demonstrated in previous 
passive resonator calculations [ 10,24] that as N ct/ is increased the round-trip losses approach a constant 
value with a periodic ripple about this value. The peaks of the ripple occur at approximately integer 
values of N ty where two different modes cross. As one of these peaks is approached, the losses of the 
current lowest-loss mode are decreasing. At the peak, both modes exhibit the same round-trip losses, and 
as the peak is passed the mode with decreasing losses becomes the new lowest-loss mode. This periodic 
dependence of the energy may help to produce a complex dependence of the MMI on Finally, N,.,, 
affects the scale of the density variation relative to the size of the mode, which increases with increasing 
N t ., r Since the density variation is dependent on the second derivative of the heating, which results from 
the Taylor series expansion ot equation (26), the scale of the density variation is expected to be reduced 
relative to the size ot the mode as N ( , t/ is increased. Together, these four factors cause the MMI to display 
a complex dependence on N eq . 

The complicated dependence of the MMI on N cq is borne out in the plots of beam quality and peak 
density variation in Figures 5 and 6, respectively. The beam quality is defined as the ratio of the far-field 


12 



Intensity 



Phase 




Figure 4. Time evolution of intensity and phase profile just before the hard-edge output mirror with N eq 
= 20 and P = 0 3 atm. The intensity is normalized to 1 at its peak value, the phase is plotted 
relative to its value at position jc = 0, and the position is in units of the width of the output 
mirror. The labels (a), (b), (c), and (d) refer to the intensity and phase profile at time 2, 8, 12, 

and 16 p.s, respectively. 
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Figure 5. Hard-edge resonator beam quality as a function of time with P — 0.3 atm and N eq = (a) 1 .0. 
(b) 2.0, (c). 3.0, (d) 4.0, and (e) 5.0. 



(a) 

(b) 

(e) 

ISJ 


0 5 10 15 20 


Time (microseconds) 


Figure 6. Peak value ot the G, profile just before the hard-edge output mirror and as a function of time 
with P = 0.3 atm and N C(/ = (a) 1 .0, (b) 2.0, (c) 3.0, (d) 4.0, and (e) 5.0. G, is in units of 
radians. This plot shows the growth of the density variation. 
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rarefaction or compression region in the medium. 

in general, the beam quality in Figure 5 remains fairly constant at about I 19 ;j» cen ^ n ^ " ^" S 

££££==! 

SSSSsSsiss 

cavity. 

Now the pressure of the gas is varied inside the cavity. Since the saturation intensity is proper- 

intend according to equation (2), causes the gain to saturate a, a tower cavity intenstty. thus reduemg 
the peak output power. 

In order to eliminate the diffraction from the edges of the output mirror, the hard-edge output 

the mode of the resonator becomes completely unstable. 

I, is noted here that variations in pressure produce results very similar to those of the hard-edge 
resonator. For the Gausstan-edge resonator, only the ^ 

edge resonator, because the relative shape of the mode does not chi an « e „' ed ear lL ' the 

number. Only the size of the mode increases with increasing Fresnel number. A * ™ e "' . F ’ , 

density variation is dependent on the second derivative of the hearing. Thus, w.lh increasing 
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Figure 7. Hard-edge resonator beam quality as a function of time with N eq = 2.0 and P = (a) 0.3, (b) 
0.5. (c) 0.7, and (d) 0.9 atm. 



Figure 8. Peak value ot the G, profile just before the hard-edge output mirror and as a function of time 
with N cq = 2.0 and P = (a) 0.3, (b) 0.5, (c) 0.7, and (d) 0.9 atm. G, is in units of radians. 
This plot shows the growth of the density variation. 




Intensity 



r ,,, r ,,, i "H i — i i rr 1 i "t-h r 




Phase 




Figure 9. Time evolution of intensity and phase profile just before the solt-edge output mirror with iV„, 
= 2.0 and P — 0.3 atm. The intensity is normalized to 1 at its peak value, the phase is plotted 
relative to its value at position _v = 0. and the position is in units of the width of the output 
mirror. The labels (a), (b), (c), and (d) refer to the intensity and phase profile at time 2. 40. 46. 
and 50 p.s. respectively. 
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numbers, the scale of the heating should be reduced relative to the size of the mode. Consequently, for 
larger Fresnel numbers a reduction is expected in the MMI effects, which is borne out in the plots of 
beam quality and peak density variation in Figures 10 and 11. With each increase in Fresnel number, the 
beam quality improves and the peak density variation grows more slowly. These trends are more obvious 
than the previous hard-edge resonator cases. 

In summary, the beam quality is much better than that of the hard-edge resonator. The initial 
value is near 100 percent compared with 19 percent for the hard-edge resonator. As mentioned before, 
the initial value is not important, since the initial beam quality can be improved by increasing the mag- 
nification [24] ot the hard-edge unstable resonator. More importantly, the beam quality remains fairly 
constant for a longer period of time than the hard-edge resonator. The beam quality begins to degrade in 
the 20- to 30-|xs range compared with the 10- to 15-|xs range of the hard-edge resonator. 


V. CONCLUSION 


A three-level lasing medium in general and COi in an unstable resonator in particular have been 
modeled. The mode-medium instability is explored numerically for various parameters in the small 
equivalent Fresnel number N,, (/ regime. The results are in agreement with those of the large Fresnel 
number regime [9], namely, the beam quality is not too sensitive to N,,, r although the larger values of N eq 
seem to provide a better beam quality for practical consideration. A lower pressure or density of the 
medium certainly improves the beam and prolongs the duration of a high quality beam at the expense of 
peak power. 

In contrast to the case of large the numerical calculations here are rather difficult, and a 
careful review of various potential problems for numerical inaccuracy are given; for example, aliasing, 
energy leakage, and Gibbs’ phenomenon. 

The major contribution of the paper is to show that the soft-edge (Gaussian reflectivity) mirror 
can indeed improve the beam quality by eliminating the hard-edge mirror diffraction ripples. The useful 
duration of the beam as well as the quality are much improved. This point is clearly demonstrated in the 
work, which also turns out to be easier to execute with better accuracy in this case. 
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f igure 10. Soft-edge resonator beam quality as a function of time with P = 0.3 atm and N eq = (a) 1 .0. 
(b) 2.0. (e) 3.0. and (d) 4.0. 



0 10 20 30 40 


Time (microseconds) 

Figure 1 1 . Peak value of the G, profile just before the soft-edge output mirror and as a function of time 
with P = 0.3 atm and N C(/ = (a) 1 .0. (b) 2.0, (c) 3.0. and (d) 4.0. G, is in units of radians. 
This plot shows the growth of the density variation. 
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